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ABSTRACT 


The  BBXGY  equations  for  M  identical,  impenetrable,  charged 
particles  which  move  in  one  dimension  and  lie  in  a  charge 
neutralizing  background,  are  shown  to  separate  into  N  uncoupled 
equations  for  the  sequence  of  N  reduced  distributions.  The 
potential  relevant  to  any  subgroup  cf  s  adjoining  particles  is 
that  of  an  s-dimensional  harmonic  oscillator  whose  frequency 
is  the  plasma  frequency  cf  the  aggregate.  The  s -particle 
spatial  equilibrium  distribution  reveals  that  particle  vibra¬ 
tions  remain  centered  about  fined,  uniformly  distributed  sites 
as  p/T  goes  from  zero  to  Infinity,  where  p  is  particle  density 
and  T  is  temperature.  Thus  it  is  concluded  that  the  system 
suffers  no  change  in  phase  for  all  p  and  T. 
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In  this  paper  we  consider  a  cne-dimensior.al  aggregate  of 
N  identical  charged  particles  which  dove  in  a  charge  neutral 
background  which  extends  ever  the  spatial  interval  (-1/2,  1/2). 
The  particles  are  further  constrained  so  that  they  are  impene¬ 
trable.  Whereas  the  equilibrium  properties  of  closely  related 
systems  have  been  studied  in  detail  [1,2],  kinetic  analyses 
have  for  the  most  part  been  restricted  to  numerical  studies  [3,4]. 
Furthermore ,  in  these  previous  investigations,  particles  were 
not  constrained  to  be  impenetrable. 

In  the  present  work,  a  general  simplifying  property  of 
the  3BXGY  sequence,  relevant  to  this  configuration,  is  derived. 
Namely,  it  is  found  that  the  hierarchy  separates  into  N  un¬ 
coupled  equations  for  the  sequence  of  N  reduced  distributions. 
The  potential  relevant  to  any  such  subgroup  of  s  labeled, 
adjoining  particles,  is  that  of  ar.  s -dimensional  harmonic 
oscillator  whose  frequencies  are  equal  to  the  plasma  frequency 
of  the  aggregate. 

The  general  solution  for  each  such  reduced  distribution 
is  obtained  and  explicit  forms  are  written  down  for  the  norm¬ 
alized  equilibrium  distributions  at  a  given  temperature.  The 
cne-particle  distribution,  for  the  first  labeled  particle,  ir. 
configuration  space,  is  found  to  be  one  half  a  gaussiar  with 
local  maximum  at  x  *  -1/2  and  local  minimum  at  x  *  +1/2.  This 
asymmetry  is  a  consequence  cf  the  constraint  of  impenetrability 
and  distir.guishability  of  particles  in  cr.e  dimension. 
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The  s-particle  equilibrium  distribution  is  found  tc  con¬ 
tain  the  Debye  distance  as  a  parameter.  Apart  from  fixed 
constanos,  the  square  of  this  length  is  equal  to  T/p,  where 
T  is  temperature  and  p  is  particle  number  density.  It  is 
found  that  as  T/p  diminishes  from  •  to  0,  the  s-particle 
spatial  distribution  passes  from  a  flat  distribution  to  one 
which  is  peaked  at  uniform  intervals.  The  exact  form  of  the 
many  -  particle  distribution  function  reveals  that  at  all  points 
of  this  variation,  particle  migrations  remain  centered  about 
sites  fixed  at  uniform  intervals.  These  findings  are  in  accord 
with  those  of  Xuns  [5],  who  studied  the  same  configuration  as 
described  herein,  through  a  partition-function  approach.  Al¬ 
though  results  are  left  in  quadrature  form,  Sunz  concludes 
that,  'the  system  is  i.;  a  crystalline  state,  at  all  tempera¬ 
tures  and  densities'.  Van  Hove  [6]  reaches  the  same  conclusion 
for  a  collection  of  finite-sized  rigid  particles  in  one  dimen¬ 
sion,  through  calculation  of  the  free  energy  of  the  configuration. 

Partition-function  analyses  of  Ler.ard  [1],ar.d  Baxter  [2] 
for  the  analogous  system  of  penetrable  charged  particles  were 
interpreted  to  imply  a  form  of  particle  binding  which  reduced 
the  effective  number  of  particles  by  one  half.  A  molecular- 
field-type  phase  transition  for  a  one- dimensional  Coulomb 
system  on  a  circle  was  reported  by  Schotte  and  Truong  [7]. 
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2.  33KGY  Equations 

We  consider  N  particles  of  charge  -  a  which  lie  in  a  one¬ 
dimensional  space  of  length  L.  The  space  is  uniformly  charged 
with  total  neutralising  charge  No.  The  Coulomb  potential 
for  this  system  [1,2]  is  given  'ey 


(2) 
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With  this  form  at  hand,  we  may  construct  the  liouviile  equation 
for  the  M-particie  distribution  (xN,p^) ,  where  x^  denotes 
the  M-di.mensional  vector,  ,  . . . and  similarly  for  ?  . 
There  results 
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Writing  di  for  dx,dp.,,  standard  reduction  [8]  of  the  second 
term  in  (3)  gives 
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The  remainder  summation  vanishes  due  to  the  boundary  conditions, 
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To  reduce  the  third  term  in  (3)  first  we  note 
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Here  we  have  written 
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Furthermore  we  are  assuming  impenetrability  of  particles  sc 
that 
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Combining  ( ^ )  and  this  last  resu 
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Chus  we  find  that  for  a  one-dime 
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fhus  the  potential  'seer.'  by  the  sub  group  of  particles 
1  through  s  is  that  of  an  s-dirensicrai  harmonic  csciliat 


with  constant  effective  spring  constant  mu  ,  varying  dis 
placement  of  or i sin,  ,  and  varying  depression  of  energy 
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The  Hamiltonian  for  this  sane  group  of  particles  is 
giver,  by  3 
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Since  che  energy  of  this  sub  group  is  separable,  the  orbits 
are  independent  [9].  Thus,  N  constants  of  the  notion  for  our 
s-particle  group  are  the  energies  (E^}.  (Here  we  are  suppress¬ 
ing  the  origin-of-tine  constant  so  that  a  single  particle  in 
one  dimension  has  only  one  constant  of  notion  [8].)  It  follows 
that  the  general  solution  to  (13)  is  given  by 

f 2 (l,...,s)  3  g(E1 , . . . ,ES )  (lo) 

where  g  is  an  arbitrary  function. 

In  equilibrium  at  the  temperature  T  one  obtains 
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**.  Asymmetry  of  Distributions 

With  f  representing  probability  density, 
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we  obtain  [together  with  the  constraint,  (3)], 
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ability  constraint  (3). 


Phase  Prooerties 


Returning  to  the  s-particle  configuration  distribution 
notion  (.21),  we  see  that  in  the  Units: 
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Note  that  the  constraint  C3 )  reappears  in  limit  (b).  To  under¬ 
stand  the  significance  of  the  Uniting  ferns  (a)  and  (b),  we 
refer  to  the  expression  for  An  given  above.  Thus,  cor¬ 

responds  either  to  T  ~  «  or  p-0,  both  of  which  extremes 
dinir.ish  the  Coulonb  interaction  between  particles  and  back¬ 
ground.  The  ir.terparticle  force  in  one  dinension  is  constant. 
The  opposite  lir.it ,  \p  -  C  ,  corresponds  either  to  T  ■*  0  or 
o-*®,  in  which  cases  the  interaction  between  particles  and 
background  is  r.ost  pronounced  resulting  in  a  rigid  spatial 
ordering  of  the  N  particles.  This  change  from  a  uniform,  fiat 
distribution  at  (p/T)  »  0  to  a  distribution  of  uniformly  spaced 
particles  at  the  limiting  value  (p/T)  *  ®,  resembles  a  contin¬ 
uous  transition  from  the  gaseous  phase  to  the  solid  phase. 
However,  the  form  (21)  indicates  that  particles  ’remember' 


their  centers  of  displacement  (x. 
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in  c/T .  In  this  sense  the  system  -air. tains  its  crystalline 
property  for  ail  o  and  T  [5], 

Conclusions 


The  33K3Y  sequence  for  N  equally  charged,  imper.etraole , 
identical  particles  which  no ve  in  a  charge  neutral  cne-dimen- 
sicr.al  space  were  shown  to  separate  into  N  uncoupled  equations 
for  the  respective  M  reduced  distributions.  The  Hamilt oniar. 
for  any  sub  group  of  connected  particles  includes  a  sumac ior.a 
cl*  ica-TCicr.ic  csci.1 cairns  vici  ccnscanc  frscusr.c^ 
"cuC  varying  di.3‘cl2.c3msn'3  cf*  E.r^c  vEiryin*  cscr*6ssicr.  c f 

energy,  lx? licit  ferns  were  obtained  for  the  coordinate  and 
ncner.tun  parts  of  the  equilibrium  s-particle  reduced  distribu- 


The  spatial  component  of  this  distribution  was  found  to 
contain  the  lecye  distance  as  a  parameter,  or  equivalently 
o/T.  In  varying  this  parameter  from  »  to  C  it  was  found  that 
the  s-cartiole  coordinate  distribution  cassed  from  a  constant 


(.flat'  function  of  s  coordinates  tc  a  product  of  s 


rune  c iens  ovsr*  cris  LangTii  or* 

the  aggregate.  However,  at  all  points  of  this  change,  particl 
coordinate  variations  remain  centered  about  what  might  be 
termed  'lattice  sites'.  Thus  one  may  conclude  that  there  is 
no  change  in  phase  of  the  system  at  any  o  or  T. 

The  key  constraint  of  this  analysis  is  that  of  impenetra¬ 
bility  tf  particles.  In  two  and  three  dimensions  £11], this 
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constraint  has  little  consequence.  However  in  one  dimension 
we  see  that  impenetrability  of  particles  implies  distinguish- 
ability  of  particles.  Thus,  reduced  distributions  are  always 
of  labeled  particles  and  kinetic  theory  becomes  a  clarifying 
tool  for  uncovering  statistical  properties  for  such  configura¬ 
tions  . 


[1]  A.  Lenard,  J.  Math.  Phys.  2,  682  (1961). 

[2]  R.  Baxter,  Proc.  Cambridge  Phil.  Soc.  59.,  779  (1963). 

[3]  J.  Dawson,  Phys.  Fluids  5,  445  (1962). 

[4]  E.  Lieb  and  D.  C.  Matt is ,  Mathematical  Physics  in  One 

Dimension,  Academic  Press,  New  York  1966. 

[5]  H.  Kunz,  Ann.  of  Physics  85,  303  (1974). 

[6]  L.  Van  Hove,  Physica  XVI,  137  (1950). 

[7]  K.  D.  Schotte  and  T.  T.  Truong,  Phys.  Rev.  A  22_,  2183  (1980) 

[8]  R.L.  Liboff,  Introduction  to  the  Theory  of  Kinetic  Equations, 

Krieger  Publishing,  Huntington,  New  York,  1979. 

[9]  H.  Goldstein,  Classical  Mechanics,  Addison  Wesley,  Cam¬ 

bridge,  Mass.  1951. 

[10]  M.  Abramowitz  and  I.  Stegun,  Handbook  of  Mathematical 

Functions,  Dover,  (Section  7),  New  York,  1965. 

[11]  R.  L.  Liboff,  J.  Math.  Phys.  7,  999  (1966). 

Acknowledgments 

This  research  was  supported  in  part  under  contract 


AFOSR- 78-3574. 


7  AuThO«C»; 

D.  M.  Heffernan  and 
R.  L.  Liboff 


9.  PERFORMING  ORGANIZATION  NAME  ANO  ADORESS 

Cornell  University 
Ithaca,  NY  14853 


I  I.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Directorate  of  Physics /W  P 
AFOSR 

Bolling  AF  Base,  D.C.  20332 


14.  MONITORING  AGENCY  NAME  &  ADDRESS^//  different  from  Controlling  Ottice) 


16.  DISTRIBUTION  STATEMENT  (o t  thin  Report) 


Approved  for  Public  Release;  Distribution  Unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  it  different  from  Report) 


19.  KEY  WOROS  < Continue  on  reverse  aide  if  necessary  end  identify  by  block  number) 

BBKGY  eqs.  One- dimens ion.  Charged  particles.  Coulomb 
interaction.  Reduced  distributions.  Plasma  frequency. 


20  ABSTRACT  /'Continue  on  reverse  aide  If  necessary  end  identify  bv  block  number) 


The  BBKGY  equations  for  N  identical,  impenetrable,  charged 
particles  which  move  in  one  dimension  and  lie  in  a  charge 
neutralizing  background,  are  shown  to  separate  into  N  uncoupled 
equations  for  the  sequence  of  N  reduced  distributions.  The 
potential  relevant  to  any  subgroup  of  s  adjoining  particles  is 


SECURITY  C  L  ASS1  F 1 C  A  Tt  Q  lyq  F  THIS  PA<3Eftt7>*n  Dm(m  Entmrmd) _ 

20.  (cont.) 

that  of  an  s- dimensional  harmonic  oscillator  whose 
frequency  is  the  plasma  frequency  of  the  aggregate. 

The  s-particle  spatial  equilibrium  distribution  reveals 
that  particle  vibrations  remain  centered  about  fixed, 
uniformly  distributed  sites  as  p/T  goes  from  zero  to 
infinity,  where  p  is  particle  density  and  T  is  temperature. 
Thus  it  is  concluded  that  the  system  suffers  no  change 
in  phase  for  all  p  and  T. 
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